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The temperature phase transition in the A'^-component scalar field the- 
ory with spontaneous symmetry breaking is investigated in the pertur- 
O"! bative approach. The second Legendre transform is used together with 

"^ ' the consideration of the gap equations in the extrema of the free energy. 

Resummations are performed on the super daisy level and beyond. The 
k> , phase transition turns out to be weakly of first order. The diagrams beyond 

^ I the super daisy ones which are calculated correspond to next-to-next-to- 

leading order in 1/A^. It is shown that these diagrams do not alter the 
phase transition qualitatively. In the limit N goes to infinity the phase 
transition becomes second order. A comparison with other approaches is 
done. 



1 Introduction 

The temperature induced phase transition in the N-component scalar Field the- 
ory (O(N)-niodel) has been investigated by many authors. The interest to this 
problem results, in particular, from the importance of the spontaneous symmetry 
breaking and symmetry restoration at finite temperature in various modern field 
theoretical models of elementary particles including different content of scalar 
fields. The present status of this problem is characterized by a general opinion 
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that the phase transition in the 0(N)-models (including N=l) is of second or- 
der (see, for instance, the standard text books |jl|, |^, |^.) This conclusion results 
mainly from the non perturbative methods - lattice simulations, average action 
(flow equation) methods, large-N expansion [|, |^, P, 1^, ||. In opposite, a first order 



phase transition was found in the most perturbative approaches [^, [T^, 0. Here, 
the key problem - the necessity of resummation of the infrared divergences of a 
series in the coupling constant - has been overcome by applying the super daisy 
resummation. However, after that another problem - the lack of the smallness of 
the effective expansion parameter near the phase transition - remains. This ob- 
servation prevents conclusive results from the resummed perturbation theory and 
it leaves room for the expectation that a further resummation of the perturbation 
expansion is able to weaken the first order character of the expansion making it a 



second order at the end. In recent papers |]T2|, |T^ by an auxiliary-mass method it 
was found that the phase transition in either the 0(1)- or the 0(A^)-model is of 
second order. Moreover, it was also stated in Ref. [|1^ that in perturbation theory 
at the two-loop level the phase transition is also of second order when diagrams 
beyond the super daisy type are taken into account. Because of the importance of 
these results we shall discuss them in more detail in what follows. Here, we note 
that the effective potential calculated by this method has an imaginary part in 
its minimum which means an inconsistency of the calculation procedure adopted. 
In our previous paper ^^ within the 0(l)-model we developed as a new 



method the combination of the gap equation resulting from the second Legendre 
transform with the condition of the free energy considered as a function of the 
condensate to be in an extremum. This allows for the exclusion of the conden- 
sate value from the gap equation and results in considerable simplifications. For 
instance, in the super daisy approximation we obtained very simple explicit for- 
mulas showing clearly a first order phase transition. We applied our method to 
a certain class of graphs beyond this approximation and obtained an increase of 
the first order character of the transition. 

In the present paper we generalize our method to the 0(A^)-model. Again, in 
the extrema of the free energy, we exclude the condensate from the gap equations 
which form a set of two equations for the Higgs and Goldstone field masses. We 
consider in detail both, the super daisy resummation and the next class of graphs 
beyond. We find that the phase transition is of first order for any finite N turning 
into a second order one in the limit of A^ goes to infinity. At the same time, the 
observed phase transition is very close to a second order one and differs from it 
in the next-to-leading order in 1/A^ for finite but large A^. We demonstrate the 
effectiveness of our method by calculating a certain class of graphs beyond the 
super daisy approximation in the 0(A^)-model, namely those which correspond 
to the next-to-next-to-leading order for large A^. Here we are left with a mostly 
numerical approach to the gap equations confirming the qualitative stability of 
our previous results. The diagram mentioned in Ref. [0 is also included in the 



carried out resummations. Its influence is discussed in the last section. 



The paper is organized as follows. In the next section we collect the basic 
formulas of the resummation method based on the second Legendre transform 
together with our method to consider the gap equations in the extrema of the 
free energy. In the third section we consider the case corresponding to the super 
daisy approximation. In the fourth section the summation of the bubble chain 
diagrams is carried out. Conclusions are drawn in the last section. Some technical 
material is banned into the appendix. 

2 Basic formulas 

We consider a scalar A^-component field theory in (3 + 1) dimensions in the 
Euclidean space-time. The action reads 

S[<P] = /dx(^i0(x)K0(a;)-^(0(x)rj, (1) 

where K = □ — rri^ is the kernel of the free action. The sign is chosen so that the 
vacuum Green functions are given by the functional representation 

Z = e^ = j D<i)e^ , (2) 

where W is the functional of the connected Green functions. It is connected by 
the relation 

W = -FT (3) 

with the free energy F. We consider the theory at finite temperature T in the 
Matsubara formalism so that the loop integrations are given by 

with the momentum p = {2ttTI, p) (/ G Z). The symbol is here a condensed no- 
tation for A^ scalar fields, = {0i, 02, . . .}. In the action Eq.(|ip the following sum- 
mation over the internal indices is assumed: 0^ = J2a=i 4>1. and 0"^ = \Yl,a=i 4>a) ■ 
In the present paper, we perform a standard renormalization procedure at zero 
temperature. 

The second Legendre transform is introduced as follows. Consider as an ar- 
bitrary functional argument the full propagator /5, instead of being the free one 

A = -K-\ (5) 

Then the functional VT is a functional depending on (3: W[l3\. Here and in the 
following we denote a functional dependence by squared brackets in opposite to 
the dependence of a function on its argument. So, for example, in momentum 



representation, /3 is a function of the momentum p: /3(p), Whereas W[P] is not, 
of course. In this sense the action Eq.(^ is also a functional of 0: S[(f)]. 
Now, the second Legendre transform results in the representation 

W = S[0] + ^Tr\ogP-^TrA-^P + W2[P] (6) 

of the connected Green functions, where f3 is subjected to the Schwinger-Dyson 
(SD) equation 

(3-\p) = A-'-J:[P]ip). (7) 

Here, S[/?](p) is the functional of all self energy graphs with no propagator inser- 
tions, 

m(P) = 2^. (8) 

where S is the functional derivative, see, e.g., |]15[ for details. PF2[/5] is the sum 
of all two particle irreducible (2PI) graphs taken out of the connected Green 
functions, 

W2[P]= E w[P] (9) 

2PI— graphs 

with the function f]{p) on the lines. 

The above formulas are written in condensed notations. So, they are valid for 
any theory with any number of fields. Below, we will write down them explicitly in 
the 0(A^)-model. But before doing that we introduce the spontaneous symmetry 
breaking by turning the sign of the mass term in the free propagator, m^ —>■ —m? 
in A so that it reads now A = p^ — m^ in the momentum representation. Then, 
on the tree level, the minimum of the energy is on the shifted fields, 0i ^ t] + v, 
where v is the field condensate. After that we change the notations according to 
{01, 02, • ■ ■} ^ {^Z + w, 01, 02, • • •} where r] is the Higgs field and {0i, 02, . . .} (a = 
1,2,. ..,A^ — 1) are the Goldstone fields which are symmetric under the residual 
0{N — 1) symmetry. 

In terms of the new variables the action reads 

/[ 7T7 A A 



-r7(n-;x^)r7 + -0(n-/i^)0 
-^ i"^' + ^''''^ + 2 (r^^ ^ 2r]v) 0^ + 0^) | (10) 

with firj = —rn? + 3-^f^ and ji^j, = —m^ + j^v'^. Again, the summation over the 
internal indices is assumed. In momentum representation, the corresponding free 
propagators are 

A^=p^ + fil and A^=/ + /i2. (11) 



The propagator for the Goldstone fields is diagonal, {A^)ab = Sab A,^. With these 
notations, the representation (j^) for the connected vacuum Green functions takes 
the form 

W = S[0] + ^Tr log P, + ^^^Tr log P^ 

-^TrA,^/3, - ^^TrA^'p^ + W2[(3,,P^l (12) 

where W2[Pn,P(i>] is the sum of all 2PI graphs with the propagators /?,, and P^j, 
on the lines corresponding to the fields t] and 0. In the 0(A^)-model, the SD- 
equations form a set of A^ equations which are equivalent to one equation due 
to the symmetry. In the considered case of the broken symmetry, two equations 
remain: 

p-\p) = A-^-S,[/?„/3^](p), 

P^\p) = A^i-E4/3„/3^](p), (13) 

where following (P) the self energy functionals are given by 

SW2[P„P^] 



^r,[P„P4{p) = 2- 



6Pr,{p) 



^ r^ ^u ^ 2 6W2[P„P^] ,^.. 



Here, the simple relation (S^)^^ = ^abS^ = 2 J^? = ^atj^^^ was taken into 
account. 

At sufficiently high temperature the symmetry is restored. This means that 
the value of the condensate is zero, v = 0, and the masses of the Higgs and the 
Goldstone fields are equal, M^ = M^. Then the SD-equations ([T3|) reduce to one 
equation, Eq. (|^. At low temperature, the symmetry is broken and v is non 
zero. The masses are different and we have to consider two equations, Eqs. ([131) . 
Then, the free energy is the function of the condensate v. It has a minimum at 
f > which can be seen already on the tree level. It may have more extrema. In 
that case due to the continuity of the free energy as a function of v an additional 
extremum for f > must be a maximum indicating a first order phase transition. 

Following ||l^ , we consider the free energy in its extrema. We take the deriva- 
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tive of W with respect to the condensate v 

'^wlR R} "^^ ^2 3A A(Ar-l) ^^ ^ dW^jP^.P^] 

In general, in this expression there are contributions proportional to dPrj/dv'^ 
and dP(i,/dv'^. But they vanish by means of the SD-equations (|T^). In Eq. (|T^, 
^^ denotes the derivative with respect to the explicit dependence on v entering 
through the vertex factors. 



Demanding now the derivative given by Eq. ([T^) to vanish 

dW2[(3r,,(3(t>] 



dv' 
we obtain the following equation, 



0, (16) 



A 2 2 ^A A(iV-l) ^^ , ^dW^ 

^^ = m - 3-Tr/3, ^^Tr/3, + 2^-^, (17) 

which has to be considered together with the SD-equations (|1^). 

In order to proceed it is necessary to make approximations. First of all, we 
put the momentum p in the SD-equations equal to zero, p = 0, that is, we turn 
to the corresponding gap equations. This is reasonable because we are interested 
in the infrared behavior and the gap in the spectrum is the decisive quantity. So 
we make the ansatz 

(3, = p' + Ml 

(3^ = p' + Ml (18) 

for the propagators and consider the SD equations (|13|) at p = 0. The assump- 
tion connected with this approximation is that the ultraviolet behavior of the 
propagator is mostly the same for the full solution and for (0), that the behav- 
ior for mediate momenta is not essential and that the important information is 
contained in the gap, i.e., in M^ and in M^. We note that in the super daisy 
approximation the corresponding contribution to S is independent of p and the 



approximation ([THD becomes an exact relation. 

The next approximation concerns the functional W2- Here, we restrict our- 
selves to graphs with at most two insertions of the trilinear vertex. In analytic 
notation we have 

W2[(3^,(3^] = W^'' + D + v^T + 0{v^). (19) 

Contributions proportional to f ^ and higher powers of f ^ are dropped. Further- 
more, we wrote down explicitely the lowest graphs, i.e., the graphs with one 
vertex, and introduced the notation D for the remaining ones. By construc- 
tion, the graphs in D do not contain f-dependent vertices and have at least two 
vertices. The graphs with one vertex are 




w^"" = ^ f Y ] +^ ( V ^ + 



3A . (o)^2 _ A A^- 1 ^(0) ^(0) _ X N^-l .^(o).2 /2m 

-1% J O AT % ^<t> A AT y^<t> > ' y'^^> 



^N^ '' ' 2 N ^ ^ AN 



where the sohd hne is the Higgs field propagator, /3^, and the broken fine is the 
Goldstone one, /3^. Taken as the approximation to W2 these graphs generate all 
super daisy graphs. The vertex factors are adduced in the Appendix, Table 2. In 
W2^ Eq. ( ^OD we introduced the notations 



.(0) 




^r = Tr/?, 



Sf = Tt(3^= (^ ; . (21) 

Now, let us rewrite the gap equations in this approximation. In the restored 
phase, we have to put v = and to take the corresponding functional derivatives. 
We obtain 






Ml = --^ + ^S(°) + A^^sf--^^. (22) 



Taking into account that W- = j^fj- holds at M^ = M^ the two equations 
become equal for symmetry reasons and we are left with one equation for the 
mass in the restored phase which we denote as Mr. 

In the broken phase, in the extremum of the free energy the condensate is given 
by Eq. (|l^. First, using the approximation (|T9|) we rewrite the gap equations 
(|T3|) as follows 

< = --^^ + I-' + |^i°^ (23) 

,A^-1 fo) ^ST 6D 

+ A^^Er-2.^--2-, (24) 

Mi = -m= + A,= + Ae(») (25) 



The equation (|T^ resulting from the condition of the extremum takes the 
form 

— = m^ - -Sf - A^^^sf + 2r. (27) 

Now, by using Eq. (p7|) we rewrite the first gap equation, Eq.(^), in the form 

-^^-(^-X|:)('»^-|^i"'-^^4"'^2r), (28) 



where we took the relation T = W- into account. Then, by means of Eq. 
we rewrite Eq. (^) as 

The second gap equation can be rewritten in a similar way. First we insert v 
from Eq. (0) in its first appearance and from Eq. (^) in its second one. We 
obtain finally 

In this way all necessary general formulas are collected. In the next sections we 
will proceed with more specific approximations. 

3 Super daisy approximation (SDA) 

The super daisy approximation (sometimes the corresponding diagrams are called 
also cactus or foam diagrams) is the result of summing up all graphs containing 
lines closed up over one vertex. It is known that this resummation is equivalent 
to take in the functional W2, Eq. (^, the contribution of W2^, (]20|), only. So, 
we have to put in the formulas given in the preceding section T = D = 0. 

First, we write down the functional W, Eq. (0), in the chosen approximation. 
To that end we rewrite the contributions like TrA~^/5^ by using the SD equations 
(|I3D and S^°) resp. Ej°) in the right handsides. Multiplying by Pnip) resp. /3^(p) 
and taking the trace, by means of Eq. (^) we have 



TrA-^/3,(p) = (SW) +c. 



where c is an irrelevant constant. Taking into account these relations we rewrite 
ly as 

m^ A 1 AT — 1 

Wiv^M,,M^) = —v' - -v' + -Trln/3, + ^— Trln/3^ - iy|° (31) 

with Wf^ given by Eq. (|20|) . Here we explictely denoted the arguments W 
dependes on. 

Now we turn to the gap equations. In the restored phase, we have v = and 
the masses are equal, M^ = M^ = Mr. There is only one equation following from 
Eqs. (H). Now it reads 

M2 = -m2 + A^^sW, (32) 



where in E^°) the propagator (3 = l/(p^ + M^) has to be inserted, similar to Eqs. 
(|1]). The corresponding value of W is W{0, M^, Mr). 

In the broken phase, we first note Eq. p9| ) which takes the form 

^"^ = ^A'.?- (33) 

The gap equations (|^), (0) turn into 

MJ = ^(4°'-Sf)- (3") 

Finally, we take the high temperature expansion for the functions entering. 
They are well known and we use them in the form given in our previous papers 



Vi(M) = — Tr n^ = + + ..., 

^ ' 2 ^ 90 24 127r ' 

^ " 12 47r +•••' ^^^^ 

where the corresponding masses Mr, M^ or M^ have to be inserted. In Eqs. (|35D , 
the dots denote contributions suppressed by higher powers of M/T. 
In this approximation the gap equations read, in the restored phase, 

m; = _„.^A(iV + 2)r- A(iV + 2)r 

12N " 8nN ' ^ ' 

and in the broken phase. 

Ml = 2m^ ^ — + 2M„ + M^— —, 

AT 

M| = ^^(M,-M^). (37) 

In the restored phase, Eq. ( |36D can be solved simply with the result 

A(A^ + 2)r 



Mr 



A(iV + 2)TV A(iV + 2)r2 

"^ + 77777 • (3^ 



TiN \\ 8nN J 12N 



We observe that it has a single solution for positive mass as long as T > T_, 

where 

I 12N , , 

m (39) 



A(A^ + 2 
9 



is the lower spinodal temperature. 

In the broken phase, the two gap equations (pTj) form a system of algebraic 
equations. While in the one component case there had been one equation of 
second order (a quadratic equation like Eq. ( ^61) in the restored phase), here the 
order is higher because there are two coupled equations. We proceed as follows. 
The second equation can be solved with respect to M^, 



AT 



{ XT y XT 



\ \Attn) 



2'kN 



M^. (40) 



This is a unique solution M^^Mjj) for all T. Note the special value M0(O) = 
independent of T. 



Now we insert M(j,{Mj^) from ( ^0|) into the first equation in formula 



Ml = ,,„._A(iV + 2)r^^ jAT A(iV-l)r 

This equation can be written down as a fourth order algebraic equation. But for 
a moment let us rewrite it formally by solving as a quadratic equation for M^ 
keeping M^{Mrj), 



3AT , 

M„ = ± . 



f^XTY , A(A^ + 2)T2 A(iV-l)T_,_, , , 



For T = we have a solution for the upper sign which is clearly the mass in the 
minimum of the free energy. Hence the solution with the lower sign describes a 
maximum. This mass vanishes for T = T^ (taking into account Eqs. (|39|) and 

Eq. (^Tj) has an explicit solution, which is, however, too large to be displayed 
here. It can be easily plotted and is shown in Figure |I|. We have chosen A = 1 
in order to make the details near the phase transition better visible. But we 
can obtain the basic properties of this solution in the follwing way. Assume we 
plot the rhs. versus Ihs. of Eq. (0). We observe that for small T there is one 
solution. At T = T_ a second solution appears and at some temperature T+ both 
solutions merge and disappear. So T+ has to be interpreted as the upper spinodal 
temperature. In this way we observe a first order phase transition. Once we have 
an algebraic solution of the gap equations, we can give an algebraic expression 
for T+. Again, it is too complicated to be displayed. However, for small A, T+ 
can be expressed as 

^ = l + t^A + 0(A2), (43) 

where the numbers tN are algebraic expressions in A^. Some numerical values are 

9 

167r- 



given in Table 1. For A^ = 1 we reobtain with ti = j^ the known result from 
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N 





1 


2 


3 


4 


5 


In 


0.0569932 


0.029511 


0.0203949 


0.0158314 


0.0130797 


0.0112318 



Table 1: Some numerical values of the coefficients t^ in Eq. 



the one component case in [Q . In the limit of large A^ we obtain 

9A 1 

tN 



167r2 (2iV) 



(44) 



which shows that the phase transition becomes second order for large A^. 




Figure 1: The masses M^, M^ and Mj. as functions of the temperature T for A = 1, 
N = 2. The upper parts of the curves for M^ and M^ show the corresponding 
masses in the minimum of the free energy, the lower parts show the masses in the 
maximum. 

Now, having solved the gap equation, we consider the functional W. Let Wb 



be its value in the minimum (broken phase), given by Eq. (^) with the masses 
inserted from the solutions discussed above and the condensate v from Eq. (|33D . 
Let Wr = W{0, Mr, Mr) be the corresponding value in the restored phase. Again, 
the analytic expressions are too large to be displayed. So, we restrict ourselves to 
some numerical values. For example, for A = 1 and A^ = 2 we have T_ = 2.44949 
and T+ = 2.52514. In T = T_ we have Wb = -8.15081 and Wr = -8.14568 
(note that the free energy differs in sign from W, see Eq. (^). In T = T+ the 
corresponding values are Wb = —8.15081 and Wr = —9.1671. In Figure |] we 
plotted F = —AW/T as function of the temperature in the region of the phase 
transition, i.e., for T_ < T < T^. The temperature T^. of the phase transition 
follows from AW = M/fe - VF^ = to be T^ = 7.767. 
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-0.0005 



-0.001 



-0.0015 



-0.002 




Figure 2: The free energy F as function of t parameterizing by means of T 
T_ + t{T^ — r_) the temperature in betweeen T_ and T+. 



4 Beyond the super daisy approximation (BSDA) 

In this section we investigate the phase transition beyond the SDA. In the chosen 
approach of representing the free energy by the 2PI functional W2 an approxi- 
mation is given by chosing an approximation for W2- While we in the preceding 



section took the 3 graphs shown in Eq. (pOD we include now all graphs consisting 
of rings (necklaces) of doubled lines. These are shown in Eq. ( [7TD in the Ap- 
pendix in generic form. In the unbroken phase where all fields have the same mass 
and there are no triple vertices, these represent all graphs where the summation 
over the internal indices is to be understood. In the broken phase we have to 
consider both kinds of fields, rj and </>, and the triple vertices. The number of the 
correponding graphs becomes larger. The motivation for this choice is that in the 
unbroken phase these are all graphs in the next-to-next-to-leading order for large 
A^. The leading order is given by the tree approximation, the next-to-leading 
order by the SDA in the preceding section. 

As the expressions in this approximation beyond SDA are quite involved we 
are forced to make further approximations. 

The aim of the present section is to find out whether the results found in 
the SDA are stable with respect in the next order. For this reason we restrict 
ourselves to the consideration of the SD equations in the broken phase, Eqs. (^) 
and (ESl), which we rewrite here in the form 



Mt 






sf + 2r 



(45) 
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where the notations 



r 

9 
f 



N ST 

N 6r 

A(iV-l)2f 5/5^(0)' 
f 6D 1 6D 



(46) 



^6 13^(0) N -161340)^ 



have been introduced. The functions 



5D 



and r 



SD 



entering here, consist 



of the graphs shown in the Appendix, Eqs. (|60|) and (ff^). The corresponding 
analytic expressions are given by Eqs. ( pB| ) and (^). Now we need to take the 
high-T approximation of these functions. For this reason we consider their generic 
form. Consider, for example, the fragment 



Tr gp^(q)aA = TYgp^(q) 



{^^?i<l) 



1 + 



3Av 



(1) 



(q) 



entering Eq. 



Note, that the external momentum is zero as we consider 



the gap equations. Now we approximate the function A, Eq. (|66D, entering 
here, by a constant, namely its value at zero momentum, A ^ A\ ^g. As A is a 
dimensionless function which does not change its sign, in the remaining integral 
the behavior of theintegrand is changed only slightly. Then we take the high-T 
approximation. The quantities a, b and c can be calculated in a standard way. 



see e.g. Si(p) in Eq. (57) in [jl4 



-SAT 1 



8nN M^ ■ ■ ■ ' 


-SATA^ + 1 1 


87riV S M^ 


-SAT 2 2 



SvriV 3Mr, + M^ 



+ ... 



+ ... 



(47) 



where the dots indicate contributions suppressed by powers of M/T. The remain- 
ing expressions correspond to graphs of the basket ball type. They have been 
calculated frequently. In the high-T approximation they read 



-SA. 



Tr,/5^(g)a = Trg/3^(g)— — Sj^)(g) 



Ti gP^{q)b = Tigprjiq)- 



N "" 
-SAA^ 



-SAT^ 

S27r2A^ 



^(1) 



N 
-2A. 



si^^(g) 



7i + • • • 
-SAT2 A^ + 1 



(1)^ 



TTMq)c = i:rMq)^^^'-;{q) 



2,2-n^N S 
SAT2 2 

74 + . . . 



72 + • • • 



S27r2A^ S 



IS 



-3A 



SAT^ 



with 



Tr,/3<^(g)a = 


= Trg/50(g)- 


AT "^ ^^) - 32vr^iV '^ ' • • • 


TrMq)b = 


= Trg/50(g)- 


^ 3 -0 WJ 327r2iV 3 


Tr,/3<^(g)c = 


= Trg/50(g)- 


-2A^(i), , -3AT2 2 , 




71 = 


. l-21n^^^ 




72 = 






73 = 


/i 




74 = 


- 1 21n2*^' + *^*, 



73 + 



(4J 



/^ 



(49) 



where /x is a normahzation constant. Again, by the dots higher orders in M/T 
are indicated. 

In this way, we have analytic expressions for the functions entering the gap 
equations (H^) which can be written as 



-3A^T 



2t^2 



327r2A^2 



{7i(l + 3A) + (Ar + l)72(e + 5) 

^,M2 + ^ + B)AB 2 ^ ^- 
a 1-AB 3^ '' 1-c 



(50) 



and 



6T 



6(3,(0) 



-A 

"aT 



:i + 2ad„ 



a(l + 3A) + 3 



{2 + a + B)AB 
l-AB 



+2 



A^ 



1 c^ 



A^ + l 1-c 






3A2T2 I (AT +1)2 



327r2A^2 j^ 3 
+ iN-l] 
+ (Ar-l) 



-73 



1 



b N -2 /I 

-^ + 2- — - - + 



(51) 



A^ + 1V3 A^ + 1-26 



1 + A 74 (2 + eA + ]B)AB 

74 + 



3 

272 c 



3ea 



1-A5 



(52) 



3 1-c. 
with a, b and c given by Eqs. (P7|), for the remaining notations see the Appendix. 
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The functions S'-^-' entering the gap equations have been given for high T in 
the preceding section by Eqs. (0). Here, in order to be consistent with (^81) , we 
include one more term into the expansion, i.e., we take 

(0) _T' MT M^L^. + 2m' 
^ "12 17+ 16^ +••• ^^^^ 

with the notation 

where 7 is Euler's constant gamma. 

Inserting these approximations we rewrite the gap equations (^5]) in the form 

f AT^ SAT \ 



XT -Ml \Ln 

27rNM -ri + M^ Sir'^Ng 



Ml = ;^ .. ""K. +^ M^ + f (55) 



with 






h 



r 



l + 2fl^ 



- _ 1 3A Ar + 2 
'^^ - ^"Svr^iV 3 ' 

The quantity D^, entering the first gap equation, depends on m^ which must be 
found from the second equation as a function of m^. 

Now, we investigate quahtatively these equations and plot both sides in one 
plot. The Ihs. is simply part of a parabola. The rhs. have a more complicated 
behavior. For large M^ resp. M^ they remain smaller than the Ihs because of 
additional factors A in front. For small M^ in the rhs. of the first equation we 
note 

Q\ 1 

so that it takes negative values for small M^. 

It is interesting to compare with the expansion of these quantities for small 
A. This corresponds to take the first graph out of the infinite sequence (0). In 
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the SD equation this means to take the graph -^ — ^^ in addition to T.^^\ Eq. 
(|2lD , in the self energy part. As an approximation, this is in between SDA and 
BSAD and does not include the effect from summing up the infinite chain of 
doubled lines. We call this perturbative approximation beyond the super daisy 
approximation (pBSDA). So, in pBSDA it holds 

h = 1 + ... 

A^ + 2 / 3A \ 3A /, 3M„ N-1^ M„ + 2M, 

5^ = 1 H In H In — 

^ 3 I AnNj 2nN 1^3 fi 




Note, that 6^ takes for sufficiently small M^ negative values in opposite to ([57|). 



The rhs. of the first equation is shown in Figure ^ in both approximations for 
a temperature near the phase transition. The curve corresponding to the pBSDA 
starts from some finite value at M^ = 0. For small temperature (T < T_ where 
T_ has to be taken in the given approximation) it starts from positive values 
and has one intersection with the Ihs. corresponding to one extremum only of 
the free energy. In opposite, for T > T_ it starts from negative values and has 
two intersections indicating the existence of a maximum. For larger T, the curve 
lowers until its two intersections with the Ihs. merge at the corresponding T+. 
This case is just shown in the figure. This is qualitatively the same behavior as 
in the SDA case. Now consider the BSDA case. Here, due to the logarithmic 
behavior in Eq. (|57D, the curve starts for any temperature from negative values 
and we have a second extremum (which must correspond to a maximum) for any 
temperature T < T+. So, formally, we observe T_ = 0. However, these negative 
values to be realized the logarithm must be large. In M^ ~ ~^^- This is the only 
qualitative change coming in from the BSDA resummation and it is of extreme 
smallness. In general, it is a quite remarkable fact that the resummation results 
in a denominator (combinations 1 — a, 1 — b, 1 — c and 1 — AB) which does not 
have a zero in the region we are interested in. Comparing the curves in Figure 
1^, we see that the quantitative difference between the different approximation is 
quite large. 

Let us now consider the second equation. It is shown in Figure §. As one can 
see, in any case in both approximations it has a solution which is similar to the 
simple solution of the second equation in the SDA case, Eq. (|^). However, for 
some value s of the parameters, there may be three solutions as seen from the 
figure. This is the substructure of the phase transition which deserves further 
investigation. 

Having clarified the general behavior of the solutions of the gap equations, one 
can start a numerical investigation. Results are shown in backing the qualitative 
similarity to the SDA. A numerical solution in the BSDA case is shown in Figure 
^. It shows qualitatively the same behavior as in the SDA case which is shown 
in the lower plot in Figure ^ for the same value of A = 0.1. We note that the 
solution in the restored phase starts in the BSAD case at considerably smaller T. 
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In fact it starts from T = but there it is exponentially small. 




Figure 3: The Ihs. and the rhs. of the first gap equation, (|55|), for different 
approximations as function of the mass M^ for A = 0.1, T = 7.7, M^ = 0.2. 




2 M 



Figure 4: The Ihs. and the rhs. of the second gap equation. 



for different 



approximations as function of the mass M^ for A = 0.1, T = 7.7, M^ = 0.2. 



5 Conclusions 



We investigated the phase transition in the 0(iV)-model by perturbative methods. 
As technical tools we used the second Legendre transform and the method of 
considering the gap equations in the extrema of the free energy. The latter 
simplifies the expressions considerably and allows for more explicit results. 

We considered two stages of approximation. We started with the super daisy 
approximation. Here, for small coupling A, which is by means of Eq. 
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equivalent to high temperature, we obtained clearly a first order phase transition. 
This generalizes the corresponding results obtained in the one component case 



found in |TJ]. The masses of the Higgs and Goldstone particles in the broken 
phase and the mass of the field in the restored phase are shown in Figure [^. The 
correponding free energy is shown in Figure ^ It should be noticed that the 
first order character of the phase transition is extremely weak. As can be seen 
from Figure |^, the change in the free energy between T^ and T_ is numerically 
small as compared to the mass scale which is set to unit in that figure. Also, the 
difference between the two spinodal temperatures, Eq. (|43|), is of order A. So, to 
leading order in A the transition is of second order and the first order character 
is a next-to-leading effect in A. For large A^, from Eq. ( ^41) it follows that the 
transition becomes second order. 

In order to check the reliability of the results obtained in SDA we calculated 
in Section 4, beyond the SDA, the next class of graphs (BSAD). In the unbroken 
phase these are all graphs in the next-to-next-to-leading order in 1/A^. They are 
given by Eq. ( [7I| ) and their dependence on A^ is taken into account completely 
(not only for large A^). These graphs can be characterised as chains of doubled 
lines. The summation of the infinite set of these graphs is performed for the 
quantities entering the gap equations. Then the resulting quantities have been 
calculated approximatively for high temperature where in essence besides the 
leading contribution only the logarithmic terms, Eq. (^Q]), have been taken into 
account. As the remaining expressions are quite comphcated, the gap equations 
have been solved numerically. An example for such a solution is shown in Figure 
1^. The general outcome is that the results obtained in SDA remain valid qualita- 
tively. They are altered in details. For instance, there exists a maximum of the 
free energy for any T < T+ so that formally we have T_ = 0. However, for small 
T this maximum is exponentially small (it is due to the logarithmic terms). Also, 
there appear additional solutions (M* in Figure |^). In this way we confirmed 
in the A^-component case the results found in |T^ for the one component case 
stating that the results in SAD are qualitatively stable with respect to higher 
loops. 

With these results obtained we are left with the difficult problem that all 
nonperturbative approaches (lattice calculations, flow equations, etc.) predict a 
second order character for the transition. In the absense of a small expansion 
parameter near the phase transition it was resonable to believe that the origin 
of the discrepancy between the perturbative and non perturbative results is due 
to this fact. If this is the case, it is important to find out the type of diagrams 
responsible for the change of the type of the phase transition. 

Now, we are going to comment some calculations where the second order phase 
transition has been determined. First we discuss the observation in Ref . [|l^ that 
the inclusion of a two-loop diagram beyond the super daisy resummation results 
in a second order phase transition. The diagram mentioned is the first one in a 
series depicted in Eq. (BUT) in the Appendix. Its exclusion does not change signifi- 



cantly the structure of the gap equations and the solutions obtained if, of course, 
all other diagrams are taken into consideration. Hence, we conclude that this 
diagram does not change the type of the phase transition. In Ref. |l7] the second 



order phase transition was observed when the renormalization group method at 
finite temperature has been applied. Results obtained in this way are difficult 
to compare with that found in the case of the standard renormalization at zero 
temperature adopted in the present paper. The point is that the renormalization 
at finite temperature replaces some resummations of a series of diagrams which 
remains unspecified basically. 

In Refs. [[T2| , |13| new non perturbative method of calculation of the effective po- 



tential at finite temperature - an auxiliary field method - has been developed and 
a second order phase transition was observed for both the one- and iV-component 
models. We left the detailed analysis of this method for other publication and 
here just mention that it seems to us not a self-consisten one because it delivers 
an imaginary part to the effective poterntial in its minima. This important point 
is crucial for any calculation scheme as a whole. Really, the minima of the effec- 
tive potential describe physical vacua of a system. An imaginary part is signaling 
either the false vacuum or the inconsistency of the calculation procedure used. 
That is well known starting from the work by Dolan and Jackiw |T^ who noted 



the necessity of resummations in order to have a real effective potential at finite 
temperature. The key point of the procedure in Refs. p!^, |13| is to relate the total 
polarization operator and the one-loop effective potential. But it is difficult to 
achive on some resonable grounds because these entities are different parts of the 
effective action in the derivative expansion. Hence, it is difficult to estimate the 
correctness of the method, its stability with respect to some small deviations, 
etc. In our method of calculations the stability of the vacuum is authomatically 
satisfied when the super daisy diagrams are resummed in the extrema of the free 
energy. 

As we have seen, the SDA possesses a lot of attractive features. In particular, 
the results obtained in this way are qualitatively stable with respect to other 
resummations. So, it is resonable to investigate in more detail its properties as 
well. 

Appendix 

In this appendix we calculate and sum up the bubble chain graphs appearing in 
the gap equations in sections 2 and 4. The corresponding well known Feyman 
rules are collected in Table 2. The vertex (0^) contains the symmetric tensor 

Vabcd = o i^abScd + 5ac5bd + ^ad^c) ■ (59) 

The functional F introduced in (|1^) consists of all graphs with 2 triple vertices, 
(?7^) or {rjcjP'), i.e., of all graphs proportional to f^ and having at least 2 vertices. 
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graphic element corresponding field propagator 



graphic element 



vertex factor notation 



6Xv 

N 



iv') 



6Xv X 



(#^ 




2A 

■ N 



iv') 



2A r 



2j,2\ 



{rfcj) 



6At/ 
— Vabcd 



N 



Table 2: The Feyman rules for the 0(A^)-model 
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In the approximation made in Section 5 the graphical representation is 



2-n 

V r 






+ ... 



+- 



1 /- 



1 

8 



1 

16 



'\ + ... 



+- 




+ 



+ 



(60) 



These are all graphs forming rings (necklaces) of doubled lines (or bubble chains) . 
We denote them by 



sS(p) 






(61) 



The loop integration in the graphs (^) are all in the just defined functions S*^^^ 
except for one integration flowing through the single line. It is useful to introduce 
the following notations, 



-3A 



-2A 



1 

N 



si^) 



(P), 



j,'<M 



(62) 



Here the vertex factors and the symmetry factor (1/2) for Til^^p) and S^ are 
included as well as the factors resulting from the summation over the internal 
indices a, 6, . . . = 1, 2, . . . , A^ — 1. For them we used 



Yl SabS, 



cd 



N 



Vi, 
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and 

/ , ^a6K6ai6iKi6ia262 • • • K„_i6„_icd^cd = (^ ^ 1) ( :^ j = Ki- 

ln terms of these abbreviations, F can be written as 

" - ~^^ TtMp) la - ^^^J ^^ TtMp) 1, (63) 



AT —i^'-'i^^' '" jY 



with 



la = - + -a^ + . . . + -eb + -eb'^ + . . . + eab + ea% + eab"^ + ... (64) 
6 2 2 2 

1 2 1 3 

76 = c2 + c3 + ..., (65) 

where the notation e = g.^'^"''-,^-, has been introduced which appears as a factor in 
front of each sequence of symbols b. 

The symmetry factors of the graphs result in factors | for sequences of sym- 
bols which are symmetric under transposition. Non symmetric sequences like ab 
appear only once because ba describes the same graph. Now we pull out the 
symmetry factor ^ and correct that for the non symmetric sequences by writing 
them twice, e.g., ab + ba instead of ab. 

In this way we are faced with the sum over all sequences of symbols a and 
b. As a symbolic notation for the subgraphs in F they have to be taken to be 
non commuting. After inserting the analytical expressions in momentum space 
representation, in fact they do commute. But first let us sum up these sequences. 

We observe that for each sequence containing a single factor a, i.e. a frag- 
ment . . . bab . . ., there is a sequence containing baab instead, baaab, etc. These 
sequences can be summed up partially, 

oo 

'^. . .ba . . . ab . . . = ... bAb . . . , 
1=1 ^^ 

where the notation 

A = 66 

1 — a 

was introduced. With other words, multiple neighbored factors a appear after this 
partial resummation only within the combination A. An exception is the sequence 
consisting of one single a which is treated separately by writing | = a — y . The 
same holds for the sequences of symbols b where in addition a factor e has to be 
taken into account. The summation results in 



"^ . . .aeb . . .ba . . . = ... aBa . 



1=1 
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with 

After these resumniations,7a, 
1 / 2a 



B 



eh 



1-6 
takes the form 



(67) 



la 



+ A + B + A^ + B' + AB + BA + . . . 



2 V 3 
Now the sequence of A's and -B's can be summed up and we arrive finally at 



_ 1 / 2a A + B + 2AB\ 
^^ ~ 2 V T ^ I- AB J 



(68) 



The sequences containing c form a simple geometric series 

cC 



oo 
i=2 



and with the notation 



we obtain 



C 



1-c 



(69) 



(70) 



7fe = c C. 

In this way, by means of Eqs. ^ and (|70D we obtained the complete expression 
for the functional F, Eq. (|63D. 

Next we have to consider the functional D appearing in W2, Eq. (p^Ol). It 
is the sum of all graphs containing only quartic vertices and at least 2 loops. 
Again, within the approximation made in Sect. 4, we restrict ourselves to the 
contributions consisting of rings of doubled lines. In the case where there is one 
field only the corresponding graphical representation is 



+ ... 





1 

48 



+ E 

n.>3 



2" n 




(71) 



which is the same as Eq. (9) in |[T^. In the A^-component case with unbroken 
symmetry the corresponding expression are essentially the same. For broken 
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symmetry they become more complicated. However, in the gap equations we 
need the derivatives SD/SPr,{p) and 5D /5[3^{p) only. By means of the relation 

5D 



6(3,{p) 



one of them is already calculated and we are left with 6D/6f3^{p). The graphs 
appearing here are quite similar to that in F, Eq. (^), where now the single line 
is P^ instead of /3^. The graphical representation of a generic graph is 



SD 



S{(3^{p))ab 




(72) 



where p is the momentum on the external line and a and b are the indices on the 
lines. Theses indices are to summed up according to 



_6D_ 
The first few graphs read then 



_5D_ 



1 

12 



N-l 
a,b=l 



6D 



5il3<pip))ab 




+ ... 




1 

+ 4 



+: 



1 

8 




1 

8 




1 
2 


1 , 



+ 



+ ... 



(73) 



We write down the corresponding analytical expression in the form 

SD 



mip) 



Ti, 



+ <?) (7i + 72 + 73) 



(74) 



+Trg/?^(p + g) 74, 
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where the 7i's correspond to the raws in Eq. ([73|). 
For 7i we have 



7i 



i^^4"MV'..i:^ 



i>2 



-6A 

N 



n+l 



^(1) 



Wi^)) Vn 



where the Vn result from the summation over the internal indices according to 



Vi = 2^ VafeaibiKii 



biab 



'n — 7 , ^abaibi ''aibia2b2 ■ ■ ■ '^a„bn 



N+l 



ab 



^)%(|)(.-., 



With the notations a, b and c given by Eq. ( |62D and b = jnijb we sum up the 
geometric series and arive at 

'b b^ 
3^T^b 



7i 



-X{N+1) 



2N 



-<--)|^r^, 



(75) 



With the remaining graphs we proceed in the same manner. For 72 we note 



1-6A /-2A^ 



72 



4 N 
+ 



N 



^(i)r^^v(i) 



^);K<i)W{<i)v, 



1 /-6A\^ /-2A\^ 



N 



N 



Tr, ( (S«(g))^ sW(g)V2 + S«(g) (sr(^))^^i) + • • • 

(76) 



-SAA^-l /I , 1 2^ 1 ,2 
iV 9 Ve e e 

which is essential the same as 7a in Eq. (|64|) except for the dependence on e. 
Summing up we obtain 



72 



-3AiV-l/ ^ i(2 + eA + i5)A5^ 

N 9 e 1 - A5 



For 73 we have simply 



73 



-A 

1^ 



{N - 1) cC 



(77) 



(78) 



so that taking all together we arrive at 



SD 



-^sTrg/3<^(p + g) 



{N + l) 
(iV-1) 
+ TvMp + q){N-l)- 



b b^ r.r ON fb b^ \ 

^ ^^^{^ + A+pA^ 
3 3 3e 1-AB 



c 



— c 



(79) 



where b 



N+l 



b. 
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Figure 5: The upper plot shows the numerical solution of the gap equations in the 
BSDA case as function of the temperature T for A = 0.1. Here M^, M^ are the 
masses in the broken phase in the minimum of the free energy (the corresponding 
masses in the maximum are not shown) , M^ is the mass in the restored phase and 
M* is the mass for the additional solutions. The lower plot shows for comparizon 
the corresponding solution in the SDA case. 
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